Bringing order through disorder: Localisation of errors in topological quantum 

memories 
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Anderson localization emerges in quantum systems when randomised parameters cause the ex- 
ponential suppression of motion. Here we consider this phenomenon in topological models and 
establish its usefulness for protecting topologically encoded quantum information. For concreteness 
we employ the toric code. It is known that in the absence of a magnetic field this can tolerate a finite 
initial density of anyonic errors, but in the presence of a field anyonic quantum walks are induced 
and the tolerable density becomes zero. However, if the disorder inherent in the code is taken into 
account, we demonstrate that the induced localization allows the topological quantum memory to 
regain a finite critical anyon density, and the memory to remain stable for arbitrarily long times. 
We anticipate that disorder inherent in any physical realisation of topological systems will help to 
strengthen the fault-tolerance of quantum memories. 

PACS numbers: 03.67.-a, 72.15.Rn, 05.30.Pr 



Introduction: Topological quantum memories are 
many-body interacting systems that can serve as error 
correcting codes These models possess degenerate 
ground state manifolds in which quantum information 
may be encoded. The size of the model and its energy gap 
then protects this information, preventing local pertur- 
bations from splitting the degeneracy and hence causing 
errors 0^3- However, the dynamic effects of perturba- 
tions when excitations are present are a serious problem 
for the stability of the memory Q , especially for non-zero 
temperature JtJ . Several promising schemes have been 
proposed [1, Q that suggest ways to combat this prob- 
lem with their own merits and drawbacks. In particular, 
in [9] it is shown that disorder can aid the stability of 
topological phases. Here we shed new light on the issue 
by showing that disorder in topological memories can in- 
duce Anderson localization [lOj . This exponentially sup- 
presses the dynamic effects of perturbations on exited 
states, and allows the memory to remain stable. 

In the definition of a topological memory, we require 
the following conditions. First, the stored information is 
encoded within a degenerate ground state of a system. 
Second, we require that the memory can be left exposed 
to some assumed noise model for arbitrarily long times 
without active monitoring or manipulations. Finally, 
measurement of the system after errors have occurred 
extracts both the (now noisy) contents of the memory 
and an error syndrome, allowing a one-off error correc- 
tion step to be performed to retrieve the original stored 
information. 

In topological memories, errors create anyonic excita- 
tions. Logical errors correspond to the propagation of 
these anyons around topologically non-trivial paths on 
the surface of the system. While the anyons are nor- 
mally static, a kinetic term emerges in the presence of 
a spurious magnetic field @. If this can act unchecked, 
even a single pair of anyons will cause a logical error after 



a time linear with the system size, L. This means that 
the memory is not resilient against any non-zero density 
of anyons initially present in the system. We demon- 
strate that randomness in the couplings of the code cause 
anyons to remain well localised in their initial positions. 
This enables the topological memory to successfully store 
quantum information for arbitrarily long times as long 
as the distribution of anyons is below a critical value. 
Hence, topological quantum memories can be made fault- 
tolerant against the dynamical effects of local perturba- 
tions. 



The toric code: The toric code is the simplest topo- 
logical memory. It is a quantum error correcting code 
whose code distance depends explicitly on the linear size 
of the system [lj|. The code can be defined on a two- 
dimensional L x L square lattice wrapped around a torus 
with a spin-1/2 particle at each edge. Stabiliser operators 
are defined on the spins around each vertex, v and plaque- 



tte, p, of the lattice, A v 
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stabilizer space, where information is stored, is composed 
of states which satisfy A v | £) = | £) and B p | £) = | £) for 
all v and p. Violations of these stabilizers on plaquettes 
or vertices are associated with quasiparticles known as 
anyons, with so-called e anyons on the vertices and m 
anyons on plaquettes 0]. These quasiparticles are hard- 
core bosons, though they also have non-trivial anyonic 
statistics with respect to each other. The effect of er- 
rors is to create and move pairs of anyons. Measurement 
of the anyon configuration provides the error syndrome, 
while error correction corresponds to a pairwise annihi- 
lation of anyons which does not form topologically non- 
trivial loops. In the limit of large system sizes, such cor- 
rection is always possible if the error probability is below 
a certain threshold. The system can therefore maintain 
a stable quantum memory even in the presence of a finite 
anyon density. The Hamiltonian of the toric code model 
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is given by 

Htc = — ^ J V A V — ^ JpBp. (1) 

v p 

Normally the case with J v = J p = J > for all v 
and p is considered. The Hamiltonian Htc has the sta- 
biliser space as its ground state and energetically pe- 
nalises anyon creation. As the e and m anyons are related 
by a duality transformation we can focus our study on 
the e anyons without loss of generality. 

Magnetic fields and quantum walks: Consider the 
following magnetic field perturbation to the toric code 
Hamiltonian, 

H h = hY J <- (2) 

i 

This perturbation commutes with the B p operators so 
it has no effect on the m anyons. However, it does not 
commute with the A v so it can create, annihilate and 
transport e anyons. It is instructive to separate this per- 
turbation into two terms = T + C with 

T = PnH h P n , C=J2 PnH h P m . (3) 

Here P n is the projector onto the subspace of states with 
n anyons of type e, satisfying J2 n Pn = 1- The operator 
T transports anyons, while the operator C is responsible 
for their creation and annihilation and hence is energet- 
ically suppressed by the gap. As such, to the first order 
of perturbation theory, the effect of Hh can be well de- 
scribed by T alone when h <C J ■ To the same order, a 
general field will also act on the e anyons in this way. 
The resulting perturbed Hamiltonian describes quantum 
walks , which transport the e anyons around the lat- 
tice. It should be noted that the perturbed Hamiltonian 
can also be mapped to the 2D Ising model in a transverse 
field, in which the spin polarized phase corresponds to 
the toric code. To proceed let us initially consider only 
the e anyons, ignoring the m's and the degeneracy as- 
sociated to them. We then deal with a reduced Hilbcrt 
space of dimension 2 L _1 , spanned by the states describ- 
ing the anyonic occupation of each of the L 2 vertices. 
The Hamiltonian of this system can be written as 

H e =^2M v ,v>t v<v > + U^2,n v {n v - 1), (4) 

v : v f v 
M v y — 5{ vv >}hi + 5 v yJ v , (5) 

where 8i v y\ = 1 only when v and v' are neighbouring 
and i is the corresponding link, and n v are the num- 
ber operators for each vertex. The operator t vv > maps a 
state with an anyon on the vertex v to one with the anyon 
moved to v' . Furthermore, it annihilates any state with- 
out an anyon initially on v. Since the anyons are hardcore 
bosons, we are interested in the case of U —> oo. 



If the initial state of the system has no anyons, the 
Hamiltonian H e has trivial effect. With no anyons to be 
moved, no logical errors can occur. Realistically, how- 
ever, single spin errors will always occur, leading to an 
ever present population of anyons. When the probability 
of such errors is low, the resulting anyon configuration 
will be of sparsely distributed anyon pairs. As such each 
pair can be considered independently, and so the study 
of two anyon walks is sufficient to determine the lifetime 
of the memory. 

The Hamiltonian for the two anyon walks is the projec- 
tion of H e onto states of two walkers only. This Hamil- 
tonian is denoted M 11 and acts upon a L 2 (L 2 — l)/2 di- 
mensional space. This is a significant reduction from the 
2 L _1 dimensional space of the full Hamiltonian, mak- 
ing numerical studies tractable. The Hamiltonian in Eq. 
(|4j) accounts only for the e anyons and ignores their non- 
trivial anyonic braiding with the m anyons that can be 
initially present in the system. This should not affect 
our results, since it has been shown that the effects of 
the statistics only arise to higher order in perturbation 
theory than we consider Q. It has also been shown in 
[l2| that, for a one dimensional walk of e anyons around 
a regular pattern of m's, the walk is equivalent to that 
when no m's are present. We would therefore expect a 
similar result to apply here, though the random pattern 
of m's will likely cause the walk to become even slower. 

In the extreme case, where only a single pair of anyons 
is present a logical error will occur when the anyons walk 
a distance of order L/2 from each other. The quan- 
tum walks induced by the field cause this to occur very 
quickly, in a time linear with the system size, L. A similar 
argument holds for an anyon density of p = N/ L 2 . Here 
traversing more than a distance of i/p/2, half the aver- 
age distance between pairs, will cause error correction to 
become ambiguous and prone to fail. Hence the dynamic 
effects of such perturbations destroy fault-tolerance when 
any non-zero density of anyons is present Q , in stark con- 
trast to the case with no field. However, by considering 
the effects of Anderson localization, fault-tolerance can 
be regained. 

Disorder and localisation: When the Hamiltonian 
is disordered, with the J not uniform but taking ran- 
dom values J v on each vertex, the theory of localisation 
places an exponential bound on the eigenstates of the 
single walker Hamiltonian [lpj . It is expected that this 
effect persists for the case of two interacting particles [3] , 
such as the anyon pairs with hardcore repulsion as con- 
sidered here. We formulate the bound on the eigenstates 
of the two particle Hamiltonian, M 11 , as follows: When 
localization takes place, each eigenstate of M 11 is peaked 
around a different position state for the anyons | vo,v' ). 
Using this fact, we denote the eigenstates | E Vo y o ) ac- 
cording to their peak. The amplitude for any other posi- 
tion state I v, v 1 ) is exponentially suppressed by the dis- 
tance of v and v' from the peak values. We define this 
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distance by, 

d(v,v';v ,v' ) = 

vam{d{v; v ) + d(v';v' ), d(v; v' ) + d(v'; v )}, (6) 

where d(v; v' ) is the Euclidean distance between v and v' , 
etc. This exponential suppression can then be expressed 



\(v,v'\E Vo y )\ < exp 



d(v,v';v a ,v' ) 



21 



vo,W 



(7) 



The width of the exponential bound on the eigenstate 
| E Vg V > o ^ is controlled by the corresponding localisation 
length, l Vo . v ' g - The localization length of the entire Hamil- 
tonian M 11 is taken to be the maximum over all these. 

From the bound on the eigenstates, a bound on 
the motion of the walkers can be derived. Taking 
Ptyfyv'fjVijV'ijt) to be the probability that the walkers 
have moved from the vertices Vi and v[ to Vf and v'f in a 
time t, we find, 

P (v f , v' f , v u v[,t) = \( Vf , v' s \e- iuUt \ Vi , v[)\ 2 



""•»o * (v f , v' f | E Vo y o ) (E Vo iV , Q | Vi, v' 4 ) I 
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V \{v f ,v' f \E„ 



\ v O,Vq 



)\\{E Vo y o \ Vi ,vi)\ 



£ eXP 

\Vo,v' 



d(wi,u-;wo,Uo) + d(vf,v' f ;v ,v' ) 
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(8) 



This sum is 

d{vi,v[;v ,v' ) 



dominated by the terms for which 
+ d(vf,v'j:;vo,v' Q ) = d, where d = 
since this is the minimum possible value. 
There are of the order of d 2 possible Vq , v' that lie directly 
on the paths between u,, v[ and Vf,v'f. Hence the bound 
can be simplified to P(v/, v'j, v t , v\, t) < d A exp[— d/l]. We 
may also note that there are of order d 3 pairs of vertices 
Vf,v'f a distance d from any initial positions Vi,v[. The 
probability P(d) that an anyon pair moves a distance d 
from its initial position therefore satisfies, 



P(d) = 0(d 7 exp[-d/Z]). 



(9) 



From this bound the exponential suppression of anyons 
moving from their initial positions is established. 

In order to verify that disorder truly does induce An- 
derson localization for the two walker Hamiltonian M 11 , 
and to determine the average value of the localization 
length for reasonable system sizes, numerical studies have 
been performed. This is done for a specific case of disor- 
der in which the J v randomly take values between J— 7/2 
and J + 7/2 according to a uniform distribution. Here J 
is the average value of the couplings and 7 is the strength 
of the disorder. To calculate the eigenstates of M 11 , the 
specific values of J and 7 need not be known. In fact it is 
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FIG. 1. The average value of the localization length I as the 
disorder y/h is increased, with each point calculated over 100 
samples. The different curves show the results for different 
system sizes L, with the units of I corresponding to the lattice 
spacing. 



only the ratio y/h, of disorder strength to magnetic field 
strength, that is important. The only constraint on J is 
that it must be large enough to maintain an energy gap 
in the full toric code Hamiltonian. It must therefore be 
at least greater than 7, but also large enough to overcome 
the reduction in the gap due to effect of the perturbation 

In our numerical study, once a Hamiltonian M is 
generated from the given disorder model, exact diagonal- 
ization is used to determine its eigenstates. To calculate 



l Vo v i for each eigenstate E Vq V i , first the peak amplitude 
and hence the values of vq and v' , are found. The prob- 
ability as a function of the distance d(v, v']Vo,v' ) is then 
derived from which l Vo y can be calculated. The localiza- 
tion length I of the full Hamiltonian M 11 is then found 
as the maximum of the l Vo y ■ We consider lattice sizes 
from L = 8 to L = 11, for which exponentially local- 
ized eigenstates can be clearly seen when y/h > 50. It is 
likely that localization also occurs at lower disorder, but 
it cannot be resolved at these lattice sizes. The resulting 
values of I for a range of disorder strengths are shown in 

Fig. m 

Localisation and error correction: The localising 
effect of the disorder is of great benefit to the toric code 
quantum memory. Suppression of anyon motion results 
in suppression of logical errors. Rather than failing in 
linear time the localisation allows the memory to remain 
stable in the presence of a finite anyon density, p, as we 
now demonstrate. Consider a regular configuration of 
anyons distributed as follows. The lattice is partitioned 
into squares of side length A, and a single anyon pair is 
created in each square. In the absence of anyon motion 
the memory remains stable, since error correction can be 
performed by simply by re-annihilating the pairs. When 
a magnetic field is applied, however, anyons can move out 
of their squares and the information of the initial pairing 
is lost. Error correction in this case can be achieved by 
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FIG. 2. The behavior and approximate values of the critical 
anyon density p c as a function of the localisation length, I. 
This is calculated by numerically solving for the side length A c 
at which the probability p of anyons moving between squares 
reaches the critical value p c where error correction fails. 



considering the parity of the number of anyons in each 
square. Initially, each square has even parity. A single 
anyon moving between neighbouring squares changes the 
parity of both. The configuration of odd parity squares 
may then be used to determine how to best undo the 
errors caused by the anyon motion. This is possible so 
long as the probability of an anyon moving over each 
boundary between neighbouring squares, p, is less than 
the critical value p c ~ 0.11 [l[. For the case of no disor- 
der, the anyons move quickly and p will exceed p c in time 
linear with A. In the presence of disorder each anyon pair 
will obey the bound of Eq. . As such, there will exist 
a square size A c for which p is kept below p c , and hence 
the memory will remain stable, for arbitrarily long times. 
This gives a critical density of p c = A^ 2 . Similarly, in the 
case that anyon pairs are created by random spin errors, 
we may expect the standard error correction procedure 
to succeed reliably when the average distance between 
pairs is at least A c , leading to the same critical density. 

It is important to determine how the critical side length 
A c and the corresponding critical density p c behave as a 
function of the localisation length, I. To do this, first note 
that the probability of an anyon moving to a neighbour- 
ing square requires d to be at least half the square size, 
A/2, giving p < J2d>\/2 Using Eq. ^ it can then 

be deduced that p = 0((X/2) 7 exp[-A/(2Z)]). The criti- 
cal side length A c is that for which p — p c . This can be 
approximated by the solution to (A c /2) 7 cxp[— A c / (21)] = 
Pc since the coefficient on the asymptotic behaviour of 
p will not have a significant effect. The critical density 
is then calculated according to p c = A~ 2 . The results 
of the numerical solution to this problem are shown in 
Fig. [5J These demonstrate the behavior p c as a function 
of I and give an estimate for its order of magnitude. It 
is found that A c = 0(1) and p c = 0(l~ 2 ), as one might 
expect from the exponential factor of Eq. ©. 



Conclusions: Magnetic field perturbations on the 
toric code induce quantum walks of anyons, which quickly 
destroy any stored information when anyons are present. 
However, we have shown that disorder induces exponen- 
tial localization which suppresses the anyon motion. This 
allows the memory to remain stable even when a finite 
anyon density is present. Since disorder will be inherent 
in any experimental realisation of topological systems, 
e.g. with Josephson junctions [14j, the effect described 
here is expected to play a significant role in their be- 
haviour. Localization will also protect against Hamilto- 
nian perturbations in other topological models, includ- 
ing those of non-Abelian anyons. The prospect of pur- 
posefully engineering disorder into topological systems to 
benefit from further localization effects, for both coherent 
and incoherent errors, is a subject of continuing study. 

Additional Note: Complementary results have been 
obtained independently by Cyril Stark, Atac Imamoglu 
and Renato Renner |15| . 
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